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ABSTRACT 


This thesis discusses the modification of and additions 
to an existing Optimal Systems Control FORTRAN Progran 
(OPTSYS) originally cktained from Professor Arthur E. Bryson 
of Stanford University. This program has been subsequently 
redesigned to run interactively on the IBM 3033 VM/CMS by 
Lieutenant Commander John G. Hoden, and additiors by 
Commander Harry A. Diel provide the user with a highly accu- 
rate graphic time response to a system designed using the 
OPTSYSX frogran. 

The addition of the FORTRAN program OPTGRAPH gives the 
user the capability tc obtain classical analysis (Pole~Zero 
Map, Root Locus, Bode, Nyquist, and Nichols) of transfer 
functions, calculated for a system designed with the OPISYSX 
program. The OPTGRAPH program uses high resolution precision 
plotting software to provide the user with highly accurate 
frequency response flcts. 

This series of frograms permits users to rapidly carry 
out simulation, analysis, and design of Optimal Systems 
Contrel problems in a totally interactive mode. 
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SYMBOLS é 


State (Ns,Ns) or Output (No,No) Weighting Matrix 

Control (Ne,Nc) Weighting Matrix 

Control Gain Matrix (Nc,Ns) 

Control (No,Nc) or Noise (No,Ng) Feedforward 

Matrix 

Open~Loop Dynamics Matrix (Ns,Ns) 

Control Distribution Matrix (Ns, Nc) 

State Disturbance Distribution Matrix (Ns,Ng) 

Measurement Scaling Matrix (No,Ns) 

Estimator Gain Matrix (Ns, Nob) 

Number of Ccntrols _ 
Number of Process Noise Sources 

Number of States 

Number of Okservations or Measurements 

White Process Noise Covariance Matrix (Ng,Ng) 
White Meas. Noise Covariance Matrix (No,No) 
Steady-State Covariance Matrix of Control (Nc,Nc) 
Control Vector (Nc,1) 

White Measurement Noise Vector (No,1), with 
zero Mean and Covariance Matrix R 

White Process Noise Vector (Ng,1), with 

zero Mean and Covariance Matrix Q 


Constant Disturbance Vector (Ng, 1) 

State Vector (Ns,1) 

Derivative cf State Vector (Ns,1) 

Estimate of State Vector (Ns, 1) 

Derivative of Estimate of State Vector (Ns,1) 
Output Vectcr (No, 1) 

Measurement Vector (No, 1) 


I. INTRODUCTION 

The purpose of this thesis is to describe the modifica- 
tion and additions tc the existing FORTRAN program (OPTSYS) 
which is used in the study and application of Oftizal 
Systems Ccntrol theory. 

The OPtimal SYstems control program was originally 
developed by Hall [Ref. 1] to support his research in 
rotary-wing aircraft control systems. Later program modifi- 
cations were made by Walker [Ref. 2] and Liu [Ref. 3] of 
Stanford University, and are designated OPTSYS 4 and OPTSYS 
5 respectively. OPTSYS modifications made by Hoden [Ref. 4] 
were primarily devoted to creating a user-friendly interac- 
tive version (OPTSYSX) of the OPTSYS 4 progran. The latest 
modifications by Diel [Ref. 5] allowed the user to save 
Matrices for subseguent runs of the OPTSYSX program, and 
formed data file sets for the time response program OPTCALC. 

The intent of this thesis work was to develop an inter- 
active frogram to plot the Pole-Zero map, Root-Locus, Bode, 
Nyquist, and Nichols responses to the open loop, closed loop 
hoise, and compensator transfer functions of a State 
Variakle Control System which has been developed using 
OPTSYS Program. Mincr modifications to the OPTSYSX Progran 
were necessary to ccrstruct a data file sets for plctting 
the Pole-Zero Map, and in calculations for the Root-Locus, 
Bode, Nyquist, and Nichols plots. 

It is assumed that the user is familiar with the ktasic 
concepts cf Control Theory and Optimal Systems Design. The 


symbol/naming conventions of Bryson [Ref. 6]. are used in 


the discussion of program operations and descriptiors of 
problems using the OFISYS Systen. 


An overview of the OPTSYSX Program capabilities and of 


modifications to the existing programs is presented first. 
This is followed by a description of the program (OPTGRAPH), 
which was developed te plot the Pole-Zero Map and to perform 
calculations for the Root-Locus, Bode, Nyguist, and Nichols 
plots. 

This work concludes with examples of various types of 
problems demonstrated in the interactive mode, including a 
copy of each terminal session with the final results. 
Complete program listings for the OPTSYS EXEC frogran, 
CPTSYSX progran, and CPTGRAPH program are included in 
appendices A,B, and C respectively. 
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A. BACKGROUND 


CETSYSX is a double-precision, interactive FORTRAN 
program employing modern control theory analysis techniques. 
Its capabilities include the calculation of the cyen-loop 
eigensystem, and the stationary closed loop system; the 
synthesis of regulators and filters; along with pcecwer spec- 
tral density, and modal distribution computations. The 
modifications introduced to the OPTSYSX program ty this 
thesis work do net affect the program's original 


capabilities. 


Be OVERVIEW 


CETSYSX is an extremely large and complex program with 
over 3000 lines of ccde. To use the program in its small 
versicn (dimensioned for a 32 xX 32 "FH, "Gt, and = "H" 
matrices), the user must extend his Virtual Machine (VM) 
memory capacity beyond the default VM memory capacity cf 720 
kilobytes to 1024 kilcbytes. 

Any significant increase in the OPTSYSX program size, 
with the the resulting requirement for additional memory 
capacity, would cause user difficulties. For this reason 
the task of obtaining the Pole-Zero maps, Root-Locus, and 
Bode, Nyguist, Nichols plots for the open loop, noise and 
compensator transfer functions was relegated to a separate 
program (OPTGRAPH). Three data files (OPTGROL DATA, OPIGRNO 
DATA, and OPTGRCM DATA), containing the open loop, noise, 
and ccmpensator transfer functions, are formed by the OFfTSYS 
program to be transfered to the OPTGRAPH prograr. The 
OPTSYS EXEC program provides an interface between the 
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OPTSYSX rprogran, the time response programs (OPTCALC and 
OPTPLCT), and the OPTGRAPH progran. 


C. OFTSYSX MODIFICATIONS 


CETSYSX program modifications consist of the addition of 
write statements, three flags, and a short routine to inter- 
pret user inputs. 

Write statements,to input data to the OPTGROL, CPTIGRNO, 
and CPIGRCM DATA files were added to the main frogran, 
subroutine INNER, sukroutine ZEROS, and subroutine POLES. 
x System information ccnsisting of number of states (Ns), 
number of controls (Nc), number of measurements (No), number 
of process noise scurces (Ng), type transfer function 
(ITFX), Markov Parameter (IE), and two flags to be used by 
the CETSYS EXEC and the OPTGRAPH programs is inputed from 
the main program. The write statements added to the subrou- 
tine INNER, and subroutine POLES input the poles for the 
open locp and noise transfer functions, and compensator 
transfer function respectively. The zeros, the numerator - 
order, the gain, and the input and output numbers for all 
three transfer functions are obtained from the subroutine 
ZEROS. 

The first of the two flags, input with the system infcr- 
mation from the main frogram to the three data files, serves 
as a marker for the CPTSYS EXEC to locate the correct data 
line for reading system inputs. The second flag, which is 
set by the routine toe interpret us.ct inputs, serves asa 
Signal for the OPTSYS EXEC to either load the OPTGRAPH 
Frogram if any one of the three transfer functions was 
calculated by the OFISYSX progran, and abort loading the 
OPTGRAPH program if none of the three transfer functions : { 
were calculated. The third flag (ITFX) is an existing ee! 
OPTSYSX frogram flag. The flag (ITFX) is passed from the 
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subroutine TF to the subroutine ZEROS to identify the calcu- 
lation cf either an open loop, noise, or compensator - 
transfer function zeros. é jacket 
The short routine added to the main program of OPTSYSX 
reads the flags ITF1, ITF2, ITF3, and IRET. The flag sent : ae 
to the CPTGROL,  OPTGRNO, or OPTGRCM DATA files is set to oe 
allow loading of transfer function data file to the CPTGRAPH re 
program if the user selects either of the transfer function 
options number two (calculate poles, residues, and zeros - 
ITF1 = 1, ITF2 = 1, or ITF3 = 1) number three (calculate 
poles, and zeros - ITF1 = 2, ITF2 =2, or, ITF3 = 2) when 
calculating the open loop, noise, or compensator transfer 
functions. The selection of any other transfer function 
option will not calculate the transfer function zeros. The 
flag IRET is read to determine if the user desires to rerun 
a problem without exiting the OPTSYSX program. Should the 
user elect to rerun a problem without exiting OPTSYSX (IRET 
= 1), the three data files are set back to the beginning 


st eneere 


(REWIND) to accept the new problem system data. : . 


D. OFTSYS EXEC MODIFICATIONS 


The OPTSYS EXEC was written by Diel [Ref. 5] to inter- 
face between the OFISYSX program andthe time resronse 
Frograms OPTCALC and CPTPLOT. The modifications to the exec 
program were to add the option to select running the 
OPTGRAPH program to the exec menues and a_ short secticn to 
interface the OPTGRAPH program with the OPTSYSX program. 

The added section to interface between the OPTSYSX and 
OPTGRAPH programs defined loader size (LDRTBLE), file defi- 
Nitions (FILEDEF), and text libraries (TXTLIB) for the 
CPTGRAPH program. The first line of the three data files is 
read to ensure that there has been at least one of the 
three transfer functions calculated by the OPTSYSX prcegram 
prior to loading OPTGRAPH progran. 
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The OFTSYS EXEC routine is also used to send Flots to 
the VERSATEC plotter from the OPTPLOT program and the 
ti . OPTGRAPH progran. . 

’ 


III. THE OPGRAPH PROGRAM 


A. PROGFAM OVERVIEW 


CPTGRAPH is an interactive FORTRAN program which uses 
the transfer function poles and zeros calculated by OPISYSX 
to calculate, tabulate, and plot the Pole-Zero Map, Pode, 
Nyquist, Nichols, and Root-Locus responses for either the 
open loop, closed loop noise, or compensator transfer 
functions. 


1. Frogram Language 


CFIGRAPH is programmed in FORTRAN following the 
conventions of the IBM System /360/370/ FORTRAN IV language. 
CPTGRAPH has been ccmpiled and run under both FORTRAN IV 
(G1) and FORTRAN H (extended) compilers on the IBM 3033. 


2. Graphics Package 


This program uses the Display Integrated Software 
System and Plotting System Software (DISSPLA) develcped and 
distributed by the Integrated Software Systems Corporation 
(ISSCC) of San Diego, California. The DISSPLA package is a 
library cf FORTRAN sukroutines using FORTRAN IV conventions. 


3. iLikrary Subroutine 


CPTGRAPH uses the International Mathematical and 
Statistical Library (IMSL) subroutine ZRPOLY to calculate 
the roots of a polynomial equation. ZRPOLY is a double 
precision FCRTRAN sukroutine to which is input a polynomial 
with real coefficients (double precision) in terms of 
decreasing powers of the polynomial. The subroutine finds 


the rects of the polynomial and returns a double precision 


complex array containing the roots. ZRPOLY is cagatle of 
accepting polynomials greater than 0 order ani less than 101 
order. 


CPTGRAPH has one main program and 24 subroutines. 
The main program and its subroutines may be divided into 
four tasic categories: 
1) File Data Input 
2) (Interactive Data Input 
3) Plot Setup and Sequencing 
4) Calculation 
A brief and general description of the program and 
its subroutines will follow in subseguent sections. 


B. GENERAL PROGRAM OPERATION 


OFTGRAPH was written to satisfy two specific objectives. 
The first is to create a program to be used as a instruc- 
tional tool for students taking controls related courses, 
and the second is as method to assist in the analysis of 
actual large order ccntrol system problems. 

To accomplish these objectives the emphasis while 
writing this program was to make the program as_ user 
oriented as possible Fy attempting to eliminate ambiguities 
and providing features which would protect the user fron 
inadvertent wrong entries. An effort was also made to mini- 
mize the Virtual Machine (VM) memory required to execute the 
program to facilitate the analysis of a large order systen. 


Three basic methods were used to protect the user 
from an abnormal program termination in the case of an inad- 
vertent wrong entry. The first was to construct the 
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subroutines to display a summery of the user inputs and 
allowing changes prior to leaving the subroutine. For the 
case where the user input is expected to fall within a given 
range, the input is examined ané@ the user is either allowed 
to continue for an input within the range or issued an error 
or warning and reasked the question for an input outside the 
range. The final method is used by the three subroutines 
which read the user input from the screen. The subroutines 
RDINT, RDREAL,and RICHAR expect an integer input, real 
humber input, and a lcgical "yes" or "no" respectively. In 
the event that a "null" line is inadvertently entered, the 
user is issued a warning and allowed another opportunity to 
enter the the correct input before abnormally terminating 
the program. The entry of two "null" lines also offers the 
user the option to to exit the OPTGRAPH program at cther 
than normal program exit points. The subroutine RDCHAR also 
examines the input and issues a warning if the input is not 
either "yes" or "no", 


2. large Order System 


The large order control problem being considered for 
analysis by the OPTSYSX and the OPTGRAPH programs is the 
X-29A aircraft longitudinal axis backup mode systen. This 
system has a (98 X 98) "F" matrix , a (2 X 98) "H" matrix, 
and a (98 X 1) "G" matrix. 

To accommodate this large order system the OPTGRAPH 
program was dimensioned to accept a system with a maximum of 
99 states (NS), 12 controls (Nc), 12 process noise sources 
(Ng), and 12 measurements (No). The Virtual Machine (VM) 
memory requirements for the program were minimized by the 
reuse of memory locations allocated for storage of calcu- 
lated data after that data has been either tabulated or 
plotted. This methcd of memory reuse has allowed the 
OPTGRAPH program to te dimensioned to accommodate a large 
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order system and still operate with less than the 1024 
KILOBYTES VM memory required to run OPTSYSX in its small 


version. - 


C. SYSTEM/MODEL DESCKIPTION 


z The system equations used by the OPTSYSX program for the 
calculation of the transfer functions are of the state vari- ’ 
able form. The system equations are: 

" system model 

xdot = (F]#x + [G]#u + [GAM ]*w (3.1) ' 
measurement equation 

z= [H]*x + (D]*u + v (3.2) ‘ 
estimator equation 

xedot = [F]*xe + [G]*#u + [K]*(z - [H]*xe) (3.3) 
open loor transfer function 

CH]*Is(I] - CFJI-1*{G] (3.4) 
closed lcop noise transfer function 

CH}*is(IlJ - [PJi-**{ Gan] (3.5) 


compensator transfer function from measurement to input 


CC)*(is{IJ-CFJ+{G J*(C )+(K J*(H JI-*) *[K] (3.6) 3 
where 
u = Control vector (Ns X 1) 4 
18 
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w = White process noise vector (Ng X 1) 

x = State vector (Ns X 1) 

xdot = Derivative of the state vector (NS X 1) 
Y = Output vector ("No X 1) 

z = measurement vector (No X 1) 


[CC] = control gain matrix (Nc,Ns) 

[0] = control feed-forward distribution matrix (Nc, Nc) 
of the control vector (u) 

{[F] = open-loorf dynamics matrix (Ns, Ns) 
(System matrix or Plant matrix) 

[G] = control distribution matrix (Ns,Nc) 

[GAM] = state disturbance distribution matrix (Ns,Ng) 
[H1 = measurement distribution matrix (No, Ns) 
{I} = identity matrix (Ns,Ns) 


D. PROGRAM ORGANIZATION 


1. Main Program 


The main program presents a menu allowing the user 
to select for analysis either the ofen loop transfer func- 
tion, the closed loop noise transfer function, the compen- 
sator transfer functicn, or exit the OPTGRAPH program. The 
data file for the selected transfer function is examined by 
the main program to ensure that the file contains the 
required transfer function data. Should the data file be 
incomplete the program gives the user a warning message 
indicating the nature of the missing data and allows the 
user te select another option. 

The main program menu: 


OFIGRAPH 


DO YOU DESIRE TO ANALYZE: 
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1. OPEN LOOF TRANSFER FUNCTION 
2. NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFER FONCTION 
4. EXIT OPTGRAPH 

ENTER OPTION NUMBER. 


Data is read from the selected file by the main 
program and formed into three arrays and seven integer vari- 
ables. The exchange of data between the main program and 
the subroutines, and between the subroutines is by the 
extensive use of labeled common statements. The data trans- 
fered through the subroutine call statements is limited 
primarily to flags, and constants. The program moves to its 
primary sequencing sukroutine (GRAPH) upon the completicn of 
file data entry. 


2. Fzrogram Sequencing 


The subroutine GRAPH first presents the user witha 
summary cf the the transfer system received from the data 
file. The user is then asked which transfer function he 
desires to to analyze using the conventions of the OPTSYSX 
program describing the transfer functions by input number 
and output number. A menu is presented which allows the 
user to select graphical and/or tabular system response 
(Pole-Zero map, Root-locus, Bode, Nyquist, or Nichols), 
select another transfer function, Or exit the OPTGRAPH 
progran. Upon completion of a graphical and/vor takular 
response cption, or change of transfer function option, the 
program returns to the GRAPH menu to allow the user to 
choose another option or exit to the main progran. 

Subroutine GRAPH data summary and menu: 


OPIGRAPH 
THE OPEN LOOP TRANSFER SYSTEM OBTAINED FROM OPTSYS CONTAINS: 


1. BOMBER OF STATES = 4 
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NUMBER OF CCNTROLS (INPUTS) = 1 
NUMBER OF MEASUREMENTS (OUTPUTS) 
MARKOV PARAMETER = 10.0**- 6 


i] 
= 


CLEAR SCREEN TO CONTINUE 


OPTGRAFH 


AN OPEN LOOP POLE-ZEKO, ROOT LOCUS, BODE, NYQUIST, AND/ 
CR NICHOLS PLOT IS DESIRED FOR: 
INPUT # ? 


CUTPOT # ? 


OPTGRAFH 


OPEN LOCP TRANSFER FUNCTION 


INPUT # = 1 
OQUTEUT #= 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR DATA FOR: 

1, POLE-ZERO MAE 

2. ROOT-LOCUS 

3. BODE 

Q. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME TRANSFER 
FUECTION) 

7. ANOTHER TYPE TRANSFER FUNCTION /EXIT OPTGRAPH 


ENTER OPTICN NUMBER. 


3. 


Graphic and/or Tabular System Response 


The selection of a graphical and/or tabular system 


response option for an Pole-Zero map, Root-Locus, Bode, open 


loop Nyquist, or ofen loop Nichols calls the subroutines 
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PZERG, RILO, BODE, NYCST, or NICHOL respectively. While the 
response of these sukroutines is different, their crgariza- 
tion is similar. 

The subroutines PZERO, RTLO, BODE, NYQST, and NICHOL 
primarily act as sequencing and plot setup subroutines. The 
interactive gquestioring of the user for graphic/takular 
response information and the system response calculation has 
been relegated to other subroutines. The subroutines PZERO, 
RTLO, BODE, NYQST, and NICHOL are divided into two major 
sections with the first being devoted to producing a graph- 
ical output and the second to producing a tabular output. 

The user is presented a subroutine menu offering the 
cption fcr graphic response, tabular response, or exiting 
the subroutine. The program moves to either the first 
section (graphic response), second section (tabular 
response) or back to the subroutine GRAPH menu depending on 
the oftion selection. The program returns to the subroutine 
menu upon completion cf the of either the graphic response, 
or takular response section, allowing the user to select 
another response opticn, or exit the subroutine. 

Features that are common to the five subroutines 
will tre described here and features that are unique to one 
of the subroutines will be described in the follcwing 
section under the aprfrepriate response heading. 

For a graphic response, the user is given tke choice 
of two printer opticns (TEK 618, or VERSATEC). With the 
selection of the TEK 618, the plot page size is defaulted to 
11 inches by 8.5 inches, and the plot will be presented on 
the TEK 618 screen. The VERSATEC option gives the user the 
option tc create up to a 21 inch by 21 inch plot (maxinun 
for the VERSATEC printer) and causes a DISSPLA METAFILE to 
be created. A scaling subroutine (PSCALE) scales the plot 
heading, iegend, etc. in proportion to the page width 
selected. The VERSATEC option offers the the advantage of a 
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high quality print from the VERSATEC printer using the 
VERSATEC printer option in the OPTSYS EXEC. The creation of 
a DISSPLA METAFILE also gives the user (after he exits 
optsys) the opportunity to use printers other than the TEK 
618 or VERSATEC. A major disadvantage with the VERSATEC 
epticn is that the user must exit the OPTGRAPH progran 
kefore he can print a graphic system response. For this 
reason, it is recommended that the system response first be 
plotted on the TEK 618 prior to selecting the VERSATEC 
option. 

For a tabular response,the user is presented a menu 
which gives him the choice of three devices (screen, 
printer, or disk) to send the tabular output. The user is 
asked if he desires to make any changes, after the tatular 
data has been calculated. If the answer is affirmative, he 
is presented a menu and allowed to make changes and rerun 
the tabular data routine. 

Examples of these features are included in the 
following sections with the interactive examples for the 
system responses. 


E. INTERACTIVE EXABFIES 


The open-loop, closed loop noise, and compensator 
transfer functions are are handled identically by the 
OPTGRAPH program, with the only difference being the program 
and graph headings identifying the transfer function tyre. 
The frequency response programs for the Bode, Nyquist, and 
Nichols responses use identical setup and tabular data 
routines. 

A gocd overall representation of the OPTGRAPH progran 
capabilities will be demonstrated by single interactive 
terminal session examples for the pole-zero, root locus, and 


bode resfonses using selected transfer functions. Tatular 
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data examples will te provided for the Pole~Zero, Root 
Locus, and Bode responses, and graphic response examples 
will be provided the Fole-Zero, koot Locus, Bode, Nyquist, 
and Nichcls responses. 


F. CHARACTERISTIC EQUATION ANALYSIS 


The subroutines PZERO and RTLO are similar in their 
basic formats. The fole and zero locations for the selected 
transfer function inrut and output numbers are read intc two 
Single dimension arrays. The extraneous zeros, calculated 
by the OPTSYSX progran, are eliminated by comparing the 
zeros to the Markov parameter sent with the system infcrnma- 
tion from the OPTSYSX progran. For a graphical response, 
the user is asked to define the plot limits in terms of 
X-coordinates (real axis) and Y-coordinates (imaginary 
axis). Points for the plot which fall outside these limits 
will be ignored. System information (tranfer function 
input number, output number, and (DC) gain) is listed in the 
graph legend. 


1. Fole-Zero Mar 


The subroutine PZ2ZERO converts the double precision 
numbers for pole and zero locations to single precision for 
compatibility with the with the Graphics Package and flots 
the poles and zeros within the plot limits defined by the 
user for the graphical response routine. For tabular data 
response the user is given system information (numerator 
order, denominator order, and transfer function (DC) gain) 
for the selected transfer function input number and outyrut 


number, and the pole and zero locations. 
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2. Eole-Zero Mar (Interactive Example) 


ss =. SS See 


The following example of a simplified open loop transfer 
function for an airplane with an autopilot in the longitu- 
dinal mode was obtained from [Ref. 7,p.315]}. 


The aircraft system is represented by: 


(DJ = [0] 


GES 


(o) 
e 
e 
(>) 
r) 
_ 
e 


[G] = 1 0. 


a. Example cf Graphic Response (Terminal Session) 


BEGIN RECORDING OF TERMINAL SESSION 
R; T=0.01/0.02 19340: 17 
EXECUTION BEGINS... 


OPTGRAPH 


DO YOU DESIRE TO ANALYZE: 
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1. OPEN LOOP TRANSFER FUNCTION 
2. NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION - e 
4. EXIT OPTGRAPH 

ENTER OPTION NUMBER. 

, 


OPTGRAPH 


THE OPEN LOOP TRANSFER SYSTEM OBTAINED FROM OPTSYS 
CONTAINS: 


1. NUMBER OF STATES = 4 
2. \NUMBER OF CONTROLS (INPUTS) = 1 
3. NUMBER OF MEASUREMENTS (OUTPUTS) = 1 


4. MARKOV PARAMETER = 10.0**- 6 . 


CLEAR SCREEN TO CONTINUE 


wre rt 


OPTGRAPH 


AN OPEN LOOP POLE-ZERO, ROOT LOCUS, BODE, NYQUIST, 
AND? OR NICHOLS PLOT IS DESIRED FOR: 
INPUT # ? 


ead 


OUTPUT # ? 


OPTGRAPH 


OPEN LOOF TRANSFER FUNCTION 


INPCT # = 17 

OUTEFOT #= 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 


1. FOLE-ZERO MAP 
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2. ROOT-LOCUS 


3.  EODE 
4. NYQUIST (POLAR PLOT) ih wd 
5. NICHOLS 
6. ANOTHER INPUT/OUTPUT COMBINATICN (SAME 
TRANSFER FUNCTION) ee 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT Oe 
CPTGRAPH 
ENTER OPTICN NUMBER. 
; 


POLE-ZERO MAP 
OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE: 
1. GRAPHICAL OUTPUT . 
2. ‘TABULAR DATA 
3. QUIT SUBROUTINE 
ENTER OPTICN NUMBER. 
4 


POLE-ZERO MAP 
OPEN LOOP TRANSFER FUNCTION . 
PLOTTER SELECTION AND PAGE SIZE oe 
(NOTE: PAGE IS 11.0 INCHES BY 8.5 INCHES WITH 
SELECTICN OF TEK618) 
WHICH PLOTTER DO YOU DESIRE: » 
1. TEK618 
2. VERSATEC 
CHOCSE OPTION 1 OR 2 


2 ° 
PAGE SIZE (MAXIMUM = 21.0 INCHES BY 21.0 ie 
INCHES) oy 
EBEIGHT = ee 
6 D 
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powcr: 


WIDTH = 


8 
PLOTTING LIMITS FOR GRAPH “e: 
X AXIS (REAL AXIS) 
X MINIMUM = ? 


-15 oe 
X MAXIMUM = ? ane 
5 
Y AXIS (IMAGINARY AXIS) 
Y MINIMUM = ? 
-5 
Y MAXIMUM = ? 
5 


PLOTTING LIMITS FOR GRAPH 
X AXIS (REAL AXIS) 


X MINIMUM = -15.00 
X MAXIMUM = 5.00 oe 
Y AXIS (IMAGINARY AXIS) oe 
Y MINIMUM =  -5.00 
Y MAXIMUM = 5.00 


DO YOU DESIRE TO MAKE ANY CHANGES ? 
TYPE "YES" OR "NO". 
n 


PLOT HEADING 


HOW MANY LINES OF HEADING DO YOU DESIRE ? 
( 3MAX ) 
ENTER NUMBER OF LINES. 
(ENTER 0" FOR NO HEADING) 


A MAXIMUM OF 32 CHARACTERS PER LINE IS 
ALLCWED 
LINE 1 IS: 

aircraft with an autopilot in 
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LINE 2 IS: 
the longitudinal mode - . 

LINE 3 IS: : = 
ref. 7; page 315 x 
PLOT HEALING 

LINE 1: AIRCRAFT WITH AN AUTOPILCT IN : 

LINE 2 ; THE LONGITUDINAL MODE e 

LINE 3 : REF. 73; PAGE 315 


DO YOU DESIRE TO MAKE ANY CHANGES ? 

TYPE "YES" OF "NO", . 
n e 
>USING A PRE-ALLOCATED DATASET FOR UNIT FT17F001. 
>USING A EFRE-ALLOCATED DATASET FOR UNIT FT18F001. 


A DISSPLA METAFILE HAS BEEN CREATED 
CLEAR SCREEN TO CONTINUE 
—~ 
POLE-ZERO MAP 
OPEN LOOP TRANSFER FUNCTION 
DO YOU DESIRE TO MAKE ANY CHANGES TO: 
1. PLOTTER / PAGE SIZE =, 
e 
2. GRAPH LIMITS 
3. HEADING 
4. NO ADDITIONAL CHANGES - PLOT POLE-ZERO 
MAP . 
5. NO CHANGES - EXIT POLE-ZERO PLOTTING 
FOUTINE . 
ENTER OPTICN NUMBER. et 
5 = 
e 
POLE-ZERO MAP a 
OPEN LOOP TRANSFER FUNCTION es 
DO YOU DESIRE: ® 
29 a 
e 
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1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE ‘ 
ENTER OPTICN NUMBER. 
3 


OPTGRAPH 


OPEN LOOP TRANSFER FUNCTION 


INPUT # = 1 

OUTPUT # = 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 


1.  FOLE-ZERO MAP 
2. ROOT-LOCUS 


3. EODE 
4. NYQUIST (POLAR PLOT) i 
5. NICHOLS oe 


6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 


CPTGRAPH aepaen 
ENTER OPTICN NUMBER. 
7 
OPTGRAPH 


DO YOU DESIRE TO ANALYZE: 
1. CPEN LOOP TRANSFER FUNCTION 
2. NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION 
4. EXIT OPTGRAPH 

ENTER OPTICN NUMBER. 

4 
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R3; T=2.80/4.67 19:345:15 
record off 
END RECORLING OF TERMINAL SESSION rea 


The preceding example follows at the end of the 
chapter as figure 3.1. 


rE. Example of Tabular Data (Terminal Session) oe 


The entering and exiting option menus for the 
tabular data terfinal sesSion are identical to the graphic 
response terminal session option menus, and have been elini- 
nated frcem the follewing terminal session. p 
BEGIN RECORDING OF TERMINAL SESSION 
R; T=0.01/0.02 19:51:03 
EXECUTICN BEGINS... Sie we 


POLE-ZERO MAP 
OPEN LCCE TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE a 
ENTER CPTICN NUMBER. . 
2 


ea] 


POLE-ZERO MAP 

OPEN LOCF TRANSFER FORCTION lata 
DO YOU DESIRE TABULAR OUTPUT TO GO TO: oe 

1. SCREEN 

2. PRINTER 

3. DISK (OPGRAFH LISTING) 
ENTER CFTICN NUMBER. 
1 
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POLE-ZERO MAP 
OPEN LOOP TRANSFER FUNCTION 
INFUT NUMBER = 1 
OUTEUT NUMBER = 1 
TRANSFER FUNCTICN (DC) GAIN = 0.1000D+01 


DENCMINATOR ORDER = 4 
LOCATIONS | 
REAL PART IMAGINARY PART 
0.0000D+00 0.0000D+00 
0.1000D+01 0.0000D+00 
-0.2000D+01 0.3464D+01 
-0.2000D+01 -0.3464D+01 


NUMERATOR ORDER = 1 


ZERO LOCATIONS 
REAL PART IMAGINARY PART 
-0.1000D+01 0.0000D+00 


DO YCU DESIRE TO CHANGE OUTPUT DEVICE? 
TYPE "YES" OR "NC". 


POLE-ZERO MAP 
OPEN LOOP TRANSFER FUNCTION 


DO YCU DESIRE: 


1. GRAPHICAL OUTPUT 


2. TABULAR DATA 
3. QUIT SUBROUTINE 


ENTER CPIICN NUMBER. 
3 


R; T=1.60/3.23 19:53:00 


record off 


Fn Sah es a 


. END RECORDING OF TERMINAL SESSION 


3. Koot-Locus t) 


The subroutine RTLO creates a closed loop. systen 
(with negative unity feed-back) from the selected transfer ae 
function input numbers and output numbers. Transfer func- coos 
tion foles and zeros are formed into two separate pclyno- ? 
Mials by the subroutine MAKPOL. The numerator polynomial 
(zeros) is multiplied by the gains, and the two polynomials 
are added to form a single polynomial. Complex roots of the 
of the pclynomial are calculated by the IMSL library sutrou- 
tine ZRPCLY. The ccmplex roots are separated into real and 
imaginary components for either plotting or tabulating. 

The plotting routine first plots the transfer func- ; 
tion pole and zero lccations then plots the root locations 
as tkey are received from the subroutine ZRPOLY. For the | 
plotting routine the gain interval specified by the user is 
divided into 2000 evenly spaced points for the subroutine > 
ZRPOLY to calculate roots. The poles, zeros, and roots are . 
plotted within the plct limits specified by the user. 

For the tabular data routine the user is given the ae 
ception tc choose between 1 and 500 points for the sukrceutine . 
ZRPOLY to calculate roots. The tabular data output provided 
for the user is system information (numerator order, denoni- 
nator order, and transfer function (DC) gain) for the 
selected transfer function input number and output nurber, e 
the pele and zero locations, and the root locus gains and 
roots 

A positive unity feedback system may be analyzed by 


specifying a negative gain range vice a positive gain range. e 
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The following interactive example is of the orfen 
l loop transfer functicn for an aircraft with an autcpfilot in 
the longitudinal mode described in Pole-Zero interactive 


secticn. 
a. Example of Graphic Response (Terminal Sessicn) 


BEGIN RECORDING OF TERMINAL SESSION 
R; T=0.01/0.02 20:28:37 
EXECUTION BEGINS... 


OPTGRAPH 


DO YCU DESIRE TO ANALYZE: 
R 1. OPEN LOOP TRANSFER FUNCTION 
2. NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION 
4. EXIT OPTGRAFH 


i ENTER CPTICN NUMBER. 
1 
iG OPTGRAFH 
THE CFEN LOCP TRANSFER SYSTEM OBTAINED FROM OPTSYS 
CONTAINS: 
’ 1. NUMBER OF STATES = 4 
Z. NUMBER OF CCNTROLS (INPUTS) = 1 
3. NUMBER CF MEASUREMENTS (OUTPUTS) = 1 
4. MARKOV PARAMETER = 10.0*#*- 6 
) CLEAR SCREEN TO CONTINUE 
OPTGRAFH 
4 AN OPEN LOOP POLE-ZERC, ROOT LOCUS, BODE, NYQUIST, 
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AND/CK NICHOLS PLOT IS DESIRED FOR: 
INPUT # ? 


CUTEUT # ? 


OPTGRAFH 


OPEN LOCF TRANSFER FUNCTION 


INPUT # = 1 

CUTFUT # = 1 
TO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. POLE-ZERO MAP 

2. FOCT-LOCUS 

3. BODE 

4. NYQUIST (POLAR PLOT) 

S$. NICHOLS 

6. ANOTHER INPOTI/OUTPUT COMBINATION (SAME 


TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 


OPTGRAPH 
ENTER CPTION NUMBER. 
2 

ROOT-LOCUS 


OPEN LCCE TRANSFER FUNCTION 
CO YOU DESIRE: 
1. GRAPHICAL OUTPUT 
2. TABULAR DATA 
3. QUIT SUBROUTINE 
ENTER CPTION NUMBER. 
1 
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ROOT-LOCUS 
CPEN LOOF TRANSFER FUNCTION 
FLOTTER SELECTION ANT PAGE SIZE 
(NOTE: PAGE IS 11.0 INCHES BY 8.5 INCHES WITH 
SELECTICN OF TEK618) 


WHICH PLOTTER DO YOU DESIRE: 
1. TEK618 
2. VERSATEC 
CHCCSE OPTION 1 CR 2 


; 2 
PAGE SIZE (MAXIMUM = 21.0 INCHES BY 21.0 
INCHES) 
HEIGHT = 
6 
WIDTH = 
8 


PLOTTING LIMITS FOR CRAPH 
X AXIS (REAL AXTS) 
X MINIMUM = ? 


fi ~15 
‘ X MAXIMUM = ? 

i) 

Y AXIS (IMAGINAEFY AXIS) 

i Y MINIMUM = ? 

-5 

Y MAXIMUM = ? 
5 


) PLOTTING LIMITS FOR GRAPH 
X AXIS (REAL AXIS) 
X MINIMUM = - 15.00 
X MAXIMUM 5.00 


Y AXIS (IMAGINARY AXIS) 


Y MINIMUM -5.00 
Y MAXIMUM 5.00 
DO YCU DESIRE TO MAKE ANY CHANGES ? 
TYPE "YES" OR "NC", 


GAIN RANGE FOR ROOT-LCCUS PLOT 
MINIMUM GAIN = ? 


MAXIMUM GAIN = ? 
500 


GAIN KFANGE FOR ROOT-LCCUS PLOT 
MINIMUM GAIN = 0.000D+00 
MAXIMUM GAIN 0.500D+03 
DO YCU DESIRE TO MAKE ANY CHANGES ? 
TYPE "YES" OR "NC", 


FLOT HEALING 
HOW MANY LINES OF HEADING DO YOU DESIRE ? 
( 3 MAX ) 
| ENTER NUMBER OF LINES. 
(ENTER "0" FOR NO READING) 


A MAXIMUM OF 32 CHARACTERS PER LINE IS 
ALLOWED 
LINE 1 IS: 
aircraft with an autopilot in 
| LINE 2 IS: 
! the lengitudinal mode 
LINE 3 IS; 
ref. 7; page 315 


| FLOT HEALING 
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LINE 1 2: AIRCRAFT WITH AN AUTOPILOT IN 
IINE 2 : THE LONGITUDINAL MODE 
LINE 3 3: REF. 7; PAGE 315 
DO YOU DESIRE TO MAKE ANY CHANGES ? 
TYPE "YES" OR "NC", 


nN 
>USING A PRE~ALLOCATED DATASET FOR UNIT FT17F001. 
>USING A PRE~ALLCCATEL DATASET FOR UNIT FT18FO001. 


A DISSPLA METAFILE HAS BEEN CREATED 
CLEAR SCREEN TO CONTINUE 


ROOT-LOCUS 
OPEN LOOP TRANSFER FUNCTION 
DO YCU DESIRE TO MAKE ANY CHANGES TO: 
1. PLOTTER / PAGE SIZE 
2. GRAPH LIMITS 
3. HEADING 
4. GAIN RANGE 
5. NO ADDITIONAL CHANGES - PLOT ROOT-LOCUS 
6. NO CHANGES - EXIT ROOT-LOCUS PLOTTING 
ROUTINE 
ENTER CPIICN NUMBER. 
6 


RCOT-LOCUS 

OPEN LOOF TRANSFER FUNCTION 
DO YCU DESIRE: 

1. GRAPHICAL OUTPUT 

2. ‘TABULAR DATA 

3. QUIT SUBROUTINE 
ENTER CPTICN NUMBER. 
3 


old 


OFTGR AFH 


OPEN LOOF TRANSFER FUNCTION 


INPUT # = 1 

CUTEFUT # = 1 
DO YCU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
CATA FOR: 


1. POLE-ZERO MAE 

2. ROOT-LOCUS 

3. BODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 
TRANSFER FUNCTION) 

7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 


OPTGRAPH 
ENTER CFIICN NUMBER. 
7 
OPTGRAFH 


DO YOU DESIRE TO ANALYZE: 
1. OPEN LOOP TRANSFER FUNCTION 
2. NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION 
4. EXIT OPTGRAPH 

ENTER CETICN NUMBER. 

4 

KR; T=12.36/15.14 20:22:54 

record off 

END RECCEDING OF TERMINAL SESSION 
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The preceding example follows at the end cf the 
chapter as figure 3.2 . 


| k. Example of tabular data (Terminal Session) 


The entering and exiting option menus for the 
tabular data terminal session are identical to the graphic 
response terminal session option menus, and have been elimi- 
nated from the fcllcwing terminal session. 

BEGIN KFECORDING OF TERMINAL SESSION 
R; T=0.017/0.02 20:42:59 


ef EXECUTICN BEGINS... 
ROCT-LOCUS 
OPEN ILOOF TRANSFER FUNCTION 
° DO YOU DESIRE: 
ie 1. GRAPHICAL OUTPUT 
- 2. TABULAR DATA 
= 3. QUIT SUBROUTINE 


ENTER CPTICN NUMBER. 
2 


ROOT-LOCUS 
OPEN LOOP TRANSFER FUNCTION 
CO YOU DESIRE TABULAR OUTPUT TO GO TO: 
1. SCREEN 
2. FRINTER 
3. DISK (OPGRAPH LISTING) 
ENTER CPTICN NUMBER. 


1 
GAIN RANGE FOR ROOT-LCCUS PLOT 
MINIMUM GAIN = ? 
1 
MAXIMUM GAIN = ? 
5 
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GAIN RANGE FOR ROOT-LCCUS PLOT 
MINIMUM GAIN = 0.1090D+01 
MAXIMUM GAIN = 0.500D+01 
LO YOU DESIRE TO MAKE ANY CHANGES ? 
TYPE "YES" OR "NC", 


n 


HOW MANY PCINTS DO YCU WANT TO TABULATE? 
(500 IS THE MAXIMUM) 
ENTER NUMBER OF POINTS. 


5 
ROOT-LOCUS 
OPEN LOCF TRANSFER FUNCTION 
INPUT NUMBER = 1 
CUTEUT NUMBER = 1 


TRANSFER FUNCTICN (DC) GAIN = 0.1000D+01 


OPEN ICCE TF POLES ANI ZEROS 


CENCMINATOR ORDER = 4 
FCLE LOCATIONS 
REAL PART IMAGINARY PART 
0.0000D+00 0.0000D+00 
Q0.1000D+01 0.0000D+00 
-0.2000D401 0.3464D+01 
-0.2000D+01 -0.3464D+01 
NUMERATOR ORDER = 1 
ZERC LOCATIONS 
REAL PART IMAGINARY PART 
-0.1000D+01 0.0000D+00 
GAIN = 0.1000L+01 
REAL PART IMAGINARY PART 
0.7074D-01 0.0000D+00 
0.8962D+00 0.0000D+00 
-0.1983D+01 0.3441D+01 
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-0.1983D+01 


GAIN = 0.1800L5+01 

REAL PART IMAGINARY PART 

0.1453D+00 0.0000D+00 

0.7948D+00 0.0090D+00 
-0.1970D+01 0. 3422D+01 
-0.1970D+01 -0.3422D+01 
GAIN = 0.2600L5+01 

REAL PART IMAGINARY PART 

0.2576D+00 0.0000D+00 

0.6552D+00 0.0000D+00 
-0.1956D+01 0.3402D+01 
-0.1956D+01 -0.3402D+01 
GAIN = 0.3400L+01 

REAL PART IMAGINARY PART 

0.4426D+00 0.1661D+00 

0.4426D+00 -0.1661D+00 
-0.1943D+01 0.33830+01 
-0.1943D+01 -0.3383D+01 


GAIN = 
REAL PART 
0 .4285D+00 
0.4285D+00 

-0.1929D+01 

-0.1929D+01 


0.4200D+01 


IMAGINARY PART 
0.3096D+00 
-0.3096D+00 
0.3363D+01 
-0.3363D+01 


DO YCU DESIRE TO MAKE ANY CHANGES ? 
TYPE "YES" OR "NC", 


ROOT-LOCUS 


OPEN ICCE TRANSFER FUNCTION 
DO YOU DESIRE TO MAKE ANY CHANGES TO: 


1. 
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OUTPUT DEVICE 
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42 


ee a Siig Bee | Pe, 
Pe ee 


oe ec ie 
By aes wee ee ear era eS 
2 7s Soe le Moe oe Se te te 
“on - et ee 


ow et a 


2. GAIN RANGE 
3. NUMBER OF POINTS CALCULATED 
4. NO ADDITIONAL CHANGES - TABULATE DATA 2 


e 
S. EXIT ROOT-LCCUS TABULAR DATA ROUTINE 
ENTER CETICN NUMBER. 
£ 
FOOT-LOCUS e 


OPEN LOOF TRANSFER FUNCTION 

DO YOU DESIRE: 
1. GRAPHICAL OUTPUT 
2. ‘TABULAR DATA e 
3. QUIT SUBROUTINE 

ENTER CPTICN NUMBER. 


3 e 
e 
R; T=1.63/3.40 20:44:54 
record off 
END RECOFDING OF TERMINAL SESSION 
. 
G. FREQUENCY RESPONSE ANALYSIS oy 
The calculations for the frequency response of the Bode, snes. 
Nyquist, and Nichols responses are identical. The only af 
difference is in the graphical presentation of the calcu- 
lated data. The subroutines BODE, NYQST, and NICHOLS call 
the subroutine FREQ to calculate the transfer function 
Magnitude and phase fcr the frequency range specified by the rf 
user. ae 
The subroutine FREQ eliminates the extraneous zeros by ae 
comparing the zeros to the Markov parameter as described in a 
the Characteristic Equation Analysis section. The frequency 7 
response magnitude and phase are calculated using double ee 
precision complex library functions following the on 
® 
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conventions of the IBM System /360/370/ FORTRAN IV language. 
The frequency response is calculated by determining the 
product of successive numerator zeros which have been 
combined with the frequency and divided by the preduct of 
successive denominatcr poles which also have been comtined 
with the frequency. The magnitude for the complex result of 
the zeros, poles, and frequency calculations is determined 
using the FORTRAN dcuble precision library functions for 
finding the magnitude of a complex argument. The phase of 
the of the complex result is determined by separating the 
real and imaginary parts of the complex results, and finding 
the phase angle using the double precision FORTRAN library 
function arctangent. The angle is converted from radians to 
degrees and is corrected for proper quadrant by checking the 
signs of the real and imaginary parts. The normal range for 
calculating phase angles by the subroutine is from +540 
degrees to -540 degrees. This range is considered adequate 
for most application, but should these limits be exceeded 
the flot will jump with a vertical line to either +360 
degrees or -360 degrees respectively before continuing with 
the plot. 

For a graphic response the user specified frequency 
range is divided into 500 equally spaced points (on a 
logrithmic scale) for which to calculate the magnitude and 
phase. The user is given the option to select between 1 ani 
500 points for tabular data. 


1. Eode Response 


The graphic response section of the subroutine BODE 
creates two plots. The first plot is the response magnitude 
in decibels versus the frequency in radians per second. The 
Magnitude from the subroutine FREQ is converted to decibels 
and plotted with its corresponding frequency. The seconi 
plot is of the phase in degrees versus frequency in radians 
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fer second. The phase from the subroutine FREQ is plctted 
with its corresponding frequency. The frequency for both 
flots is plotted on a logrithmic scale. The scale range for 
both is set automatically to include the entire range of 
data to ke plotted by the subroutine MINMAX which -ans the 
data for the maximum and minimum values. System information 
(tranfer function ingut number, output number, and (fC) 
gain) is listed in the graph legend. 

The magnitude in decibels and the phase in degrees 
is searched for "0" decibel and 180 degree crossover feints 
to calculate the phase and gain margins. 

The tabular response section gives the user the 
opticn to select the the number of frequency points to be 
calculated between 1 and 500 for the frequency range he has 
specified. For takular data response the user is given 
system information (numerator order, denominator order, and 
transfer function (DC) gain) for the selected transfer func- 
tion input number and output number, and the frequency 


(cadians per second), magnitude (decibels), phase (degrees), 


and the real and imaginary frequency response parts. 
2. Bode 


The fcllcwing example of compensator transfer function with 


a filter and regulator was obtained from [Ref. 8,pp. 382 - 
384 j. 
The compensated system is represented by: 


[C] = [ -254.1 -19.57 ] 


[0° 
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[GAM] =10. { 
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[K} = 1 95.44 
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a. Example cf Graphic Response (Terminal Session) 


BEGIN RECCRDING OF TERMINAL SESSION — 
R; T=0.01/0.02 18:30:23 
graphics 

EXECUTION BEGINS... 


OPTGRAPH 


DO YOU DESIRE TO ANALYZE: 
1. OPEN LOOP TRANSFER FUNCTION 
2. NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION os 
4, EXIT OPTGRAPH cape 
ENTER OPTION NUMBER. : 
3 


OPTGRAPH 
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THE COMPENSATOR TRANSFER SYSTEM OBTAINED FRCM 
OPTSYS CCNTAINS: 


{ 1. NUMBER OF STATES = 2 ! 
2. NUMBER OF CONTROLS (INPUTS) = 
3. NUMBER OF MEASUREMENTS (OUTPUTS) = 1 
: 4. MARKOV PARAMETER = 10.0**- 6 : 
| CLEAR SCREEN TO CONTINUE 
OPTGRAPH 
2 t 
é 
AN OPEN LCCP POLE-ZERO, ROOT LOCUS, BODE, NYQUIST, 
AND/OR NICHOLS PLOT IS DESIRED FOR: 
INPIT # ? 
» . : 
CUTPUT # ? 
1 
- , 
‘ OPTGRAPH ee 
COMPENSATCR TRANSFER FUNCTION 
INPUT # = 1 — 
" OUTPUT # = 1 ro 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR . 
DATA FOR: 
1. EOLE-ZERO MAP 
p 2. ROOT-LOCUS 
3. EODE 


4. NYQUIST (POLAR PLOT) 
5. NICHOLS 
® 6. ANOTHER INPUT/OUTPUT COMBINATION (SAME ' 
TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 
CPTGRAPH 
e ENTER OPTION NUMBER. , 
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BODE PLOT 

COMPENSATCR TRANSFER FUNCTION ' 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 

2. TABULAR DATA 

3. QUIT SUBROUTINE ' 
ENTER OPTICN NUMBER. 
1 


BODE PLOT 
COMPENSATCR TRANSFER FUNCTION 
PLOTTER SELECTION AND PAGE SIZE 
(NOTE: PAGE IS 11.0 INCHES BY 8.5 INCHES WITH 
SELECTICN OF TEK618) 
WHICH PLOTTER DO YOU DESIRE: 
1. TEK618 
2. VERSATEC 
CHOCSE OPTION 1 OR 2 


2 
PAGE SIZE (MAXIMUM = 21.0 INCHES BY 21.0 
INCHES) 
REIGHT = 
5.65 
WIDTH = 
8.5 


PLOTTING LIMITS FOR GRAPH 
FREQUENCY RANGE - RADIANS PER SECOND 
MINIMUM FREQUENCY = ? 
0.1 
MAXIMUM FREQUENCY = ? 
1000000 


PLOTTING IIMITS FOR GRAPH 
FREQUENCY RANGE - RADIANS PER SECONT 
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MINIMUM FREQUENCY 

MAXIMUM FREQUENCY 

DO YOU DESIRE TO MAKE ANY CHANGES ? 
TYPE "YES" OR "NO", 


0. 100E+00 
0. 100E+97 


n 


PLOT HEADING 
HOW MANY LINES OF HEADING DO YOU DESIRE ? e 
( 3 MAX ) 
ENTER NUMBER OF LINES. 
( ENTER "O" FOR NO HEADING ) 


A MAXIMUM OF 32 CHARACTERS PER LINE IS 
ALLCWED 
LINE 1 IS: 
filter sinulation e 
LINE 2 IS: . 
ref. 8; pages 332 - 334 


PLOT HEADING ea 
LINE 1: FILTER SIMULATION 
LINE 2 : REF. 8; PAGES 332 - 334 

DO YOU DESIRE TO MAKE ANY CHANGES ? 
TYPE "YES" OR "NO", ° 


A DISSPLA METAFILE HAS BEEN CREATED 
CLEAR SCREEN TO CONTINUE 


GAIN MARGIN AND PHASE MARGIN 
PHASE CROSSOVER FREQUENCY IS 75.1184 RAD/SEC; e 
GAIN MARGIN IS -59.8906 DB 
GAIN CROSSOVER FREQUENCY IS ******eee*e FAD/SEC; 
FHASE MARGIN IS -89.9480 DEG 
CLEAR SCREEN TO CONTINUE e 
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BODE PLOT 
COMPENSATCR TRANSFER FONCTION 
DO YOU DESIRE TO MAKE ANY CHANGES TO: 
1. FELOTTER / PAGE SIZE 
2. GRAPH LIMITS (FREQUENCY RANGE) 
3. HEADING 
4. NO ADDITIONAL CHANGES - PLOT BODE FLCTS 
5. NO CHANGES - EXIT BODE SUBROUTINE 
ENTER OPTICN NUMBER. 
5 


BODE PLOT 
COMPENSATCR TRANSFER FUNCTION 
DO YOU DESIRE: 

1. GRAPHICAL OUTPUT 
2. TABULAR DATA 
3. QUIT SUBROUTINE 
ENTER OPTICN NOMBER. 
2 


OPTGRAPH 


COMPENSATCR TRANSFER FUNCTION 

INPUT # = 1 

OUTPUT # 1 
DO YOU DESIRE GRAPHICAL RESPONSE AND/OR TABULAR 
DATA FOR: 

1. FOLE-ZERO MAP 

2. ROOT-LOCUS 

3. EODE 

4. NYQUIST (POLAR PLOT) 

5. NICHOLS 

6. ANOTHER INPUT/OUTPUT COMBINATION (SAME 


it 
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TRANSFER FUNCTION) 
7. ANOTHER TYPE TRANSFER FUNCTION /EXIT 


CPTGRAPH 
ENTER OPTION NUMBER. 
7 
OPTGRAPH 


DO YOU DESIRE TO ANALYZE: 
1. OPEN LOOP TRANSFER FUNCTION 
2. NOISE TRANSFER FUNCTION 
3. COMPENSATOR TRANSFER FUNCTION 
4. EXIT OPTGRAPH 

ENTER OPTION NUMBER. 

4 

R; T=7.01/8.84 18:39:11 

record off 

END RECORLING OF TERMINAL SESSION 


The preceding example follows at the end of the 


chapter as figure 3.3 and figure 3.4. 
bk. Example of Tabular Data Listing File 


The terminal session for the Bode tabular data response is 
Similar to the Bode graphic response terminal session. 
Appendix D contains an example of tabular output sent toa 
disk (OPGRAPH LISTING) for five points between 10 and 100 
radians/second. 


3. Nyguist Respense 


The graphic response section of the subroutine NYCST 
creates a polar plot from the magnitude and phase data 
received from the sukroutine FREQ. The plot is automati- 
cally scaled from a minimum radius of one by using the 
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subroutine MINMAX to determine the maximum phase and magni- 
tude values. The gain / phase rargin calculations and the 
tabular data output are identical to the subroutine BODE. 


4. Nyquist Respcnse (Interactive Example) 


The following example is of the open loop transfer 
function for the filter simulation example described in the 
Bode interactive section. 

This example follows at the end of the chapter as 
figure 3.5. 


5. 


The graphic response section of the subroutine 
NICHCL creates a single rectangular plot of magnitude in 
decibels versus phase in degrees. The magnitude from the 
subroutine FREQ is converted to decibels and plotted with 
the corresponding magnitude. The plot is automaticaily 
scaled using the subroutine MINMAX to determine the maxinum 
phase and magnitude values. The gain / phase margin calcu- 
lations and the tabular data output are identical to the 
Subroutine BODE. 


6. Nichols Response (Interactive Example) 


ee 


The following example is of the open loop transfer 
function for the filter simulation example described in the 
Bode interactive section. 

This example follows at the end of the chapter as 
figure 3.6. 
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IV. CONCL 


A. CCNCLIUSIONS 


An evaluation of the computational ability of OPTSYSX 
and CPTGRAPH, the frogram was tested using an 82 X 82 
matrix, provided by NASA-Edwards, of longitudinal eguations 
for the X-29A experimental forward-swept wing fighter 
aircraft prototype. The OPTSYSX program array had to be 
redimensioned and a 2-megabyte virtual machine size was 
required for a system of this size. 

The open loop eigenvalues calculated by the OPTSYSX 
program compared favorably with the Eigensystem data 
supplied by NASA-Edwards, but the OPTSYSX program failed in 
calculating the open loop transfer function zeros. 
Fxtraneous zeros (calculated by OPTSYSX) were of the same 
approximate magnitude as the open loop system zeros. This 
prevented the identification and elimination of the extra- 
neous zeros by comparison with the Markov parameter. The 
open loop zero identification was further complicated by not 
having available open loop transfer function zero informa- 
tion from NASA-Edwards to compare the OPTSYSX results to. 

When requested for additional information NASA-Edwards 
Frovided a computer tape with a revised 98 X 98 matrix of 
longitudinal motion equations for the X-29A fighter 
aircraft. For this revised system, the OPTSYSX program 
requires a virtual machine size in excess of 2-megatyte. 
Time constraints have prevented running the OPTSYSX program 
with the revised systen. 

The OPTGRAPH program was tested using the OPTSYSX 
Frogram's open loop Eigenvalue and transfer function cutput 
data for the X-29A 82 X 82 system of longitudinal motion 
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equations. OPTGRAPH provided excellent quality Bode glcts 
for all of the oper loop transfer function inputs and 
outputs and did not exhibit any computational difficulties 
as indicated by the lack of error messages at the end of the 
test. While the test results are encouraging as to the 
ability of the OPTGRAFH program to assist in the analysis of 
large order systems, the actual computational accuracy of 
the OFPTGRAPH program can not be fully verified until accu- 
rate transfer functicn information can be obtained from the 
OPTSYSK program. 


Be. RECCMMENDATIONS 


Eased on the results of this Thesis the following areas 


emerged as possible areas for further research and study. 
1. 


The OPTSYSX program has two similar root finding 
subroutines (HQR and HQR2), to find the transfer function 
Eigensystem (poles) and zeros. The order of the transfer 
function numerator is calculated separately in the sutrou- 
tine ZEROS. 

The inclusicn of the open loop transfer function 
eigensystem (pole locations), zero locations, transfer func- 


tion numerator orders, and transfer function gains, for the 


revised (98 X 98) Ilcngitudinal mode system by NASA-Edwards 


should simplify the lccating the problem in calculating the 
transfer function zeros. 

The open loof system data from NASA-Edwards also 
included Bode magnitude and phase plots, and gain and fhase 
margin information, which can be used to evaluate more 
conclusively the OPTGRAPH program's ability to help evaluate 
a large crder systen. 
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The OPTSYSX program requires in excess of 2-megabyte 
of virtual machine memory when configured for large matrix 
operations (98 X 98). Virtual machines with this memory 
capacity are not normally available to a user. The memory 
usage for matrix storage and manipulaticn is a possible area 
for the reduction in program memory requirement size. The 
OPTSYSX program retains most of its computational arrays in 
memory. This method simplifies programming but is extremely 


inefficient in its use of virtual machine memory. 


3. Further Modifications 


eS Se ae Se eee 


The area of modern digital controls should be further inves- 
tigated. The computational abilities of the OPTGRAPH 
program make it readily adaptable, with minor modifications, 
for assisting in the analysis of transfer functions in the 
digital domains. 
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APPENDIX A 
OPTSYS EXEC LISTING 


TRACE OFF 
PES ESSER ESET EEC SE SSE ST SSE SS SS EST STS TSE SL SSE SSS EL ESS SSS SES SESS 
* THE OPTSYS EXEC 
* CONTROLS THE OPTSYSX, OPTCALC AND OPTPLOT 
* TO DETERMINE THE TIME RESPONSE OF A SYSTE”. 
* BY H. A. DIEL 
* VERSION 1.0 16 JULY 1984 
* MCDIFIED BY M. H. LAPTAS 
PE SSE RES EEE SS EES TEES SES PETS T SET SSS STL ETS SS SSS TSE TSS SS ES ESET SS 
CLRSCEN 
EBEGTYPE -ENDFOUR 


We He ee Re He He He He tee He ie he ee a ee Re a KK KEK KKK KKK KEKE ERE KE 


CUTPUT FROM OPTSYSX WILL COME TO THE SCREEN IF YOU 
WISH LTISK FILE 
ENTER OD ANY OTHER INPUT YIELD SCREEN 


WR HH He He He He He ee HB He He ee ee KK ee KK KK EK 


~ENDFCUR 

SREAD VAES 6ANS 

SIF .&ANS NE .D &GOTC -ONE 

FILEDEF 06 DISK OUTPUTX LISTING A1 
CLRSCEN 

EBEGTYPE -ENDSIX 
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PELE SESE SE SES TESTES ESTES SCTE SES SEE EES STE ST EE SET ES TS SE SES Et SS 


OUTPUT WILL GO TO DISK FILE *OUTPUTX LISTING Ai! 


RHR ERREKEK KEKE KKK EEEREKEEKKEREKE EEE EEE KKEEKEEEEKEEKSEK 


-ENDSIX 
CP SLEEP 3 SEC 

-ONE 

CLRSCRN 

CP LINK 0039P 191 251 RR 

ACC 281 F 

WEEK KKK ERA KK HE EEK EE EERE KKK KEKE KAEKKEEKEKEKKK KEE EEKEK EEK KK 
* * 
* ALICW THE USE CF AN OLD "OPTMAT DATA Ai" * 
* * 


REE KERR EKEKE KEKE BERGE ERE REE AEKKEK KEE KEREKEREK KKK EEKK EEK EK 
RENAME OPTMAT DATA A1 OPTSYS DATA A1 

SIF &RC NE 0 &GOTO -FIRST 

RENAME OFTSYS DATA Al OPTMAT DATA A1 

FINDSTAK OPTMAT DATA A1CO1 0 LIMO002 ALL GROUP1 

GREAD VAFS &A1 &A2 

SIF .&A2 NE .1 &GOTO -START 

CLRSCEN 

EBEGTYPE -ENDTWO 


HERKEEREREEAKEK AREER DEKE ERAREREKEKEKERKAEEEKKE KEEEK EEK EK 


YOU HAVE A DATA FILE NAMED "OPTMAT DATA* ON YOUR A 
CISK THAT WAS PREVIOUSLY GENERATED BY THE OPTSYS 
FROGRAM AND CCNTAINS THE F, G, H, GAMMA, A AND B 
MATRICES FROM THAT RUN. 
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IF YOU WOULD LIKE TO USE THESE SAME MATRICES FOR 
THIS RUN, THE CPTSYS PROGRAM WILL READ IN THE 
DESIRED DATA AT THE APPROPRIATE TIME, 


IF YOU TYPE (Y) ES. 
ANY OTHER INPUT KILL RESULT IN THAT FILE BEING ERASED! 


HTK RRR A RE PAK EE REE REE ERR RRA REE EK EE EK 
-ENDTWO 

EREAD VARS &ANS 

SIF .&ANS EQ .Y SGOTO ~START 


“FIRST 
MRK RHR RK RR RK RK KK KK RRR EE RK ER IK OK RK HK Ke 


ERASE THE OLD "“CPTMAT DATA Ai" DATA FILE 
FLACE “000 O" IN THE NEW “OPTMAT DATA FILE" 
TO ACT AS A FLAG FOR OPTSYSX AND OPTCALC 


+ ££ + & & 
+ + & © 


FI SOR ITO I RR RIO AIO RIO ERR RR tok ak 
ERASE OPTMAT DATA A1 

&STACK 000 0 

FILESTICK OPTIMAT DATA Al F 80 1 

-START 

CLRSCEN 

EBEGTYPE -ENDONE 


THE OPTSYS EXEC CONTROLS FOUR PROGRAMS: 


1 OPTSYSX FCRTRAN (SYSTEM ANALYSIS) 
GENERATES OPTMAT AND OPTGRAPH DATA SETS 
- RERUN USES OPTMAT 

2 OPTCALC FCRTRAN (CALCULATE TIME RESPONSE) 
GENERATES OPTPLOT DATA SET USES OPTMAT FROM 
OPTSYSX 

3 OPTGRAPH FORTRAN (POLE-ZERO, ROOT-LOCUS, ECDE, 
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NYQUIST, ETC) 
USES OPIGROL -OPTIGRNO -OPTGRCM FROM OPTSYSX 
4 EXIT 


OPTPLOT FORTRAN IS THE FOURTH PROGRAM 
USES OFTPLOT DATA SET FROM OPTCALC 
TEE SIZE OF THE DATA SETS VARY WITH THE SYSTEM 
CRDER, AND CAN USE ABOUT 20% OF THE USERS DISK SPACE. 
THEREFORE ENSURE THAT SUFFICIENT DISK SPACE IS 
AVAILABLE. 
IF DATA SET IS ALREADY AVAILABLE YOU MAY KUN ANY 
PROGRAM 
~ ENTER 1, ANY OTHER INPUT 
RETURNS TO MENU 
-~ENDONE 
SREAL VAES &ANS 
SIF .&ANS EQ .1 &GOTO -OPTSYS 
SIF .&ANS EQ .2 &GOTO -OPTCALC1 
SIF .&ANS EQ .3 &GOTO -OPGRAPH1 
SIF .6ANS EQ .4 EGEXIT &RC 
&GOTO -START 
~OPTSYS 
FILEDEF 8 DISK OPTPLCT DATA Al (PERM 
FILEDEF 9 DISK OPTMAT DATA Al (PERM 
FILEDEF 10 DISK OPTGRCL DATA Ail (PERM 
FILEDEF 01 DISK OPTGRNO DATA Al (PERM 
FILEDEF 04 DISK OPTGRCM DATA Al (PERM 
GLOBAL TXTLIB VFORTLIE IMSLDP NONIMSL 
STYPE LCADING OPTSYS.... GENERAL STATE VARIABLE ANALYSIS 
PROGRAM 
ETYPE wee. OPTMAT DATA Al ...-. INPUT DATA SET ON RERUNS 
OPTSYSX 
CLRSCEN 
STYPE..- DATA SETS .. OPTMAT AND «+ OPGRAPH CREATED 
CP SLEEP 5 SEC 


&EGOTC -START 

-OPTCALC1 

PES Se PEPE TS SPE LES SSS SSE SE SSS SLES SSS EST SS LS SS TS TS TSS SSE SE SS ST 
* * 
* CHECK FOR DATA IN THE FILE "OPTMAT DATA " * 
* BEFORE LOADING OPTCALC * 
x * 
We He Re he Re he Re ee ee he he eRe he hee ee he He eo oe a eo ee ok eo ok Ke ok KK RK 
FINDSTAK OPTMAT DATA A1C0O1 0 LIMO002 ALL GROUP1 

G&READ VARS &A1 &A2 

SIF .&A2 EQ .1 &GOTO -OPTCALC 

ETYPE PROPER DATA FILE IS NOT AVAILABLE FOR OPTCALC 

CP SLEEP 5 SEC 

&GOTO -START 

~OPTCALC 

ETYPE...ee2-e--OPTCALC.... TIME RESPONSE PROGRAM 

FILEDEF 8 DISK OPTPLCTI DATA Al (PERM 

FILEDEF 9 DISK OPTMAT DATA A1 (PERM 

GLOBAL TXTLIB VFORTLIE IMSLDP NONIMSL 

CPTCAIC 

ETYPE «+ OPTPLOT LATA A1 CREATED 

CP SLEEP 5 SEC 

CLRSCEN 

EBEGTYPE -ENDNINE 


PEEP eR ESE ESE EER EEE EEE ELE TES SESE SES TEE TEST TSS EST ST SETS ET SS SS TF 


iF YOU ARE CISSATISFIED WITH THE RESULTS 
THUS FAR AND WOULD LIKE TO EXIT TO CMS, 


-TYPE "Y¥* TO EXIT- 


(ANY OTHER INPUT TO CONTINUE) 
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LN 


PES SSS SESE ESTES SE TS SEL SE SSL SEES SES ST SESE TT SS SS TS SSS EST SE TSS 


-~ENDNINE 

SREAD VARS &ANS 

SIF .&ANS EC .¥Y &GOTC -START 

ETYPE FLOTS NEXT 

CP SIEEP 4 SEC 

GLOBAL TXTLIB DISLIBVS 92DISLIB 92INTLIB VFORTLIB GRFLIE 
NCNIMSL CMSLIB 

FILEDEF 8 DISK OPTPLOT DATA Al (PERM 

FILEDEF 11 DISK LBLANK TMP A3 (RECFM F LRECL 2400 BLKSIZE 
2400 BLKSIZE 2400 XTENT 600 

FILEDEF 12 DISK ISSCCMAP MAPDTA * (RECFM F LRECL 400 
BLKSIZE 400 

FILEDEF 13 DISK TABLET TMP A3 (RECFM VS LRECL 2068 BLKSIZE 
208 XTENT 1000 

FILEDEF 14 DISK LBLANK MAPDTA §SLBMODE (RECFM VS LRECL 608 
BLKSIZE 608 

FILEDEF 17 DISK DISSFIA SYSUT1 A3 (RECFM FB LRECL 2000 
BLCCK 2000 XTENT 500 

FILEDEF 18 DISK DISSFLA METAFILE AY (RECFM VBS LRECL 19065 
BLCCK 19069 


OPTPLOT 

EGOTC -META 

-OPGRAPH1 

Pics PRE EEE EEE TTR SS TS SER ESE LE LS SESS TT ESTES ES SE SES SEES ES SS SS SS 
* * 
+ CHECK FOR DATA IN THE OPTGRAPH DATA SETS " * 
* EEFORE LOADING OPGRAPH * 
* * 


HAH RK RH KR RR RR ER a 
HR RE AIO I RO ROK Bk eR OR OR EO A 
RENAME OCETGROL DATA A1 OPTSYS DATA Al 

éIF EFC NE 0 &GOTO -CFETNOISE 

RENAME OPTSYS DATA A1 OPTGROL DATA Al 
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FINDSTAK OPTGROL DATA A1 CO1 0 LIMO002 ALL GROuP1 

GSREAD VARS &A1 &SA2 SA3 SA4 SAS 

SIF .&§A2 NE .0 &GOTO -SIXTEEN 

-OPTNCISE 

RENAME QGETGRNO DATA A1 OPTSYS DATA Al 

SIF &RC NE 0 &GOTO ~CETCOMP 

RENAME OPTSYS DATA Al OPTGRNO DATA Al 

FINDSTAK OPTGRNO DATA A1 CO1 0 LIMO02 ALL GROUP1 

SREAD VAES &A1 6A2 &A3 &A4 EAS 

SIF .&A2 NE .0 &GOTO -SIXTEEN 

-OPTCCMP 

RENAME OPTGRCM DATA A1 OPTSYS DATA A1 

SIF &RC KE 0 &GOTO -NCDATA 

RENAME OPTSYS DATA A1 OPTGRCM DATA A1 

FINDSTAK OPTGRCM DATA A1 CO1 0 LIMO002 ALL GROUP1 

SREAD VARS &A1 &A2 &A3 SA4 EAS 

SIF .6A2 NE .0 &GOTO -SIXTEEN 

-NODATA 

EPRINT PROPER DATA FILE IS NOT AVAILABLE FOR OPGRAPH 

CP SLEEP 5 SEC 

&GOTO -START 

-SIXTEEN 

6TYPE ..eOPTGRAPH.... CLASSICAL ANALYSIS OF OPTSYS OUTPUT 

SET LDRTELS 10 

GLOBAL TXTLIB DISLIBVS 92DISLIB 92INTLIB VFORTLIB GRFLI8 
NONIMSL CMSLIB 

SET LDRTBLS 10 

FILEDEF 01 DISK OPTGRNO DATA At (PERM 

FILEDEF 02 DISK OPGRAEFH LISTING Al 

FILEDEF 03 PRINTER (RECFM FA BIKSIZE 133 PERM 

FILEDEF 04 DISK OPTGEFCM DATA Al (PERM 

FILEDEF 06 TERM (RECFM FA BLKSIZE 133 

FILEDEF 05 TERM (RECFM FA BLKSIZE 80 

FILEDEF 10 DISK OPTGEKCL DATA A1 (PERM 

FILEDEF 11 DISK LBLANK TMP A3 (RECFM F LRECL 2400 BLKSIZE 
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2400 XTENT 600 
FILEDEF 12 DISK ISSCCMAP MAPDTA * (RECFM F LRECL 400 
BLKSIZE 400 = 


69 


FILEDEF 13 DISK TABLET TMP A3 (RECFM VS LRECL 208 BLKSIZE ss 
208 XTENT 1000 
FILEDEF 14 DISK LBLANK MAPDTA &LBMODE (RECFM VS LRECL 608 
BLKSIZE 608 e 
FILEDEF 17 DISK DISSFLA SYSUT1 A3 (RECFM FB LRECL 2000 
BLCCK 2000 XTENT 500 
FILEDEF 18 DISK DISSFIA METAFILE AY (RECFM VBS LRECL 19065 
BLCCK 19069 « 
CPTGRAPH 
“META 
Pee ee ELE SES ELE SS SELES ESE SESS SESS SCT SSS SE SST SCS SES SSC SSS SESS SSL Ss 
* * e 
* CHECK FOR FILE "DISSPLA METAFILE A" ON * 
* THE USER'S DISK BEFORE GOING TO DISSPOP * 
* * 
PIF ESE ESE ST FES SELES ESTE SELES TE SETS SSS STE SST SSS TS TSS SSS TT SES SF e eas 
RENAME DISSPLA METAFILE A4 OPTSYS METAFILE A4 
SIF &RETCOLE NE 0 &GCTO -START 
RENAME OPISYS METAFILE A4 DISSPLA METAFILE A4 pe 
-EIGHT 7 
CLRSCEN 
SBEGTYPE -ENDTEN 8 
° 
e 
RAK MK REE RHE EKK KE KKK KK EERE KEK AEE KK EK KKK KEKE KEKEKKAKKAKEK ESE KK EE os 
DC YOU WANT A VRSTEC PLOTTER SMOOTH COPY OF THE ed 
THE DISSPLA METAFILE THAT YOU JUST CREATED? ® 


Penk ee, Ta! Sit ae Fs eS ee aS fae Ps Sa Ga Tce Re See, Saat Mage. ee a la an Oo gs Rag ES pt a Se Mila, See, Soo Mad hg gear Sree ema ae Cal So Bei te at pcre en tele ie ge 
sae MAT ets, Saha 8 ee acy eI Me ge es cate SSO gene ae ae ta gg WELL meee Gane, TSM eae Ee ce. ES Rye Ie aa A Tian Sart NOS, OS ea ise Sates Gg eee em oe, Siege OS, SAPS aE Sg ke cig Pings, URN 
ee gE eg SOT ge Oe Se fe Na ale NO we aT eg gi Met OT ca eT Sg ye ge Se gr gS eg ek SN eg Pe ig ge 8 gee Se Bee ae ena ere ee ee a em 


(Y OR N) 


PEELE ESE SES ESE ST EES STS SSE LE LESTE SETS ETS TS TSS TES STE TS SST TTS ET SS 


-ENDTEN 
SREAD VARS &ANS 

GIF .SANS EQ .¥ &GOTC -NINTH 
SIF .SANS EQ .N &GOTC -START - 
&GOTO -EIGHT 

-NINTH 

EXEC DISSPOP VRSTEC 

CLRSCEN 

SBEGTYPE -ENDTWELVE 


MR RR RHR RR RRR KARE RRR KR ERE ER RRR RICK FORK doe 
YOUR GRAPH(S) CAN BE PICKED UP AT THE COMPUTER CENTER. 
THE GRAPH(S) WIII BE ADDRESSED TO "POP (USER ID)". 

MAAK RR BR RR RK FRR RRR RIOR RIOR RII $a RR tk ok 


~ENDTWELVE 


CP SIEFEP 5 SEC 
&GOTO -START 
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aati 


a fe Tet 


APPENDIX B 
OPISYSX PROGRAM LISTING 


This portion of the thesis contains the OPTSYSX FORTRAN 
program (88 pages). 
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